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Abstract. Quasicrystals can carry, in addition to the classical phonon 
displacement field, a phason displacement field, which requires a generalized 
theory of elasticity. In this paper, the third-order strain invariants (including 
phason strain) of icosahedral quasicrystals are determined. They are connected 
with 20 independent third-order elastic constants. By means of non-linear 
elasticity, phason strains with icosahedral irreducible F^-symmetry can be 
obtained by phonon stress, which is impossible in linear elasticity. 



PACS numbers: 61. 44. Br Quasicrystals, 62.20.Dc Elasticity, elastic constants 

1. Introduction 

Apart from the ordinary phonon degrees of freedom, quasicrystals (QC) have phason 
degrees of freedom, referring to relative shifts of the constituent density waves ^ . 
Therefore, and due to their lack of translational order, quasicrystalline structures are 
usually constructed as irrational cut of a decorated hyperspace structure by physical 
space E\\ |, |. The phason de grees of freedom are connected with the displacement 
field along the orthogonal space E-^. The generalized elasticity is described in terms 
of spatially varying phonon and phason displacement fields |l|, |^ . 

Icosahedral QC have three phonon and three phason degrees of freedom and 
associated components of a displacement field. Linear elastic theory provides five 
independent second-order elastic constants, two belonging to pure phonon elasticity, 
two to pure phason elasticity and one to a coupling between phonons and phasons. 

Within linear phonon elasticity, icosahedral QC behave essentially like isotropic 
media |^]. Faithful icosahedral symmetry exists for physical properties described by 
tensors of rank > 5 only Accordingly, the non- linear elasticity of icosahedral 
QC is anisotropic 0. 

Fundamental research on classical non-linear elasticity has been performed many 
years ago 1^, ^ . The authors of l^j, |ll], |3| have already determined the four 
linearly independent, icosahedral elastic tensors of rank six, related to third-order 
phonon elastic invariants. Contrary to this, in the isotropic case one has only three 
independent third-order phonon elastic invariants, or elastic constants. Ishii has 
calculated the pure phason third-order icosahedral invariants. The aim of this paper 
is to generalize the classical non-linear elasticity and to determine all third-order 
icosahedral elastic invariants, which occur when phason strains are included. The 
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idea leading to this work was to find a possibility to generate phason F'^-strain in 
icosahedral QC by phonon stress, which is impossible within linear elasticity. 

2. Generalized elastic theory of icosahedral QC 

According to their icosahedral diffraction pattern, the mass density of icosahedral QC 
is a sum of density waves indexed by a reciprocal lattice L of icosahedral symmetry: 
pix) = J2keL '^^viik ■ x) = J2keL \Pk\exp[i{k ■ x + 4>k)]- The phases 0k of the basis 
vectors of L are six degrees of freedom |^] , parametrized by the phonon and phason 
displacement fields u and w via 0fc — (pkfi — fc" ■ u — ■ w . Here, fc" — k and fc-'- 
are the projections of reciprocal six-dimensional hyperlattice vectors onto i?" and , 
respectively. The phonon and phason displacement fields u and w are the projections 
of the hyperspace displacement field uQ) w onto i?" and E-^, respectively. 

We denote the position of a point in the undistorted QC a and the corresponding 
position in the distorted structure x, where x = u + a. In the Lagrangian scheme, 
all quantities depend on the variable a The phonon strain tensor rj" has its 
components of the classical symmetric form 

1 / dui duj \ 1 duk duk 



2 V daj doi I 2 dai da^ 
or, written in terms of the Jacobian Fij = 77^- = Fkj — Sij) . This strain 

tensor is free of rigid rotations. The phason displacement gradient ^ splits into a 
and a F^ part (see (^), and both are assumed to increase the elastic energy [l5|] . 
Therefore, we have a phason strain tensor 77™ with 

In the linear limit |ff^| ^ 1, the components of 77" take their well-known shape 

'lij ~ 2 V da-j ^ dai>- 

The isothermal Helmholtz free energy F{ri",ri^) per undistorted volume can be 
expanded into the Taylor series 

^ ~ 2 ^ijkl Vij Vkl + g ^ijklmn Vij Vkl Vmn ' ' ' ^ 2 ^"^^^ ~^ Q ^ij'^ Vk + ■ ■ ■ ^ (3) 

with C°'ji.i being second-order and C°'j^ij^n third-order cartesian elastic constants due 
to Brugger p6| , which is perhaps the most familiar notation (i, j, k, l,m,n G {1, 2, 3}; 
b G {u,w}). In the right part of (||), the irreducible strain components of Appendix 



a 



^ are used. Here we have, e.g., i G {1, . . . , 6} if a = u and i G {1, . . . , 9} if a = w. 
Because of the index permutation symmetries and the symmetries of the QC, not all 
of these C's are independent. If the elastic energy is to be written with independent 
elastic constants Cf** and Cf'"^ of second and third order only, one has to use the 
invariants Zf'' and If '"^ of the generalized elasticity: 

F = Cflf + Cf^lf" + ... . (4) 

These invariants (and also the expansions (|^)) must fulfill the condition I{gri"- , grj^) = 
I{ri^,rf") for any transformation g of the icosahedral group Y (or Yh). Clearly, 
symmetries like J,""™ = '"^ = ... and C™"" = C;'"'" = ... are assumed to hold. The 
elastic constants of (^ follow from appropriate repeated differentiation with respect 
to components of rf^ and rf" . 
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A generalization of the classical result shows that the generalized Piola- 
KirchhofF stresses and t™, which are measured in the undistorted state, have 
components 

^i] - Q^a — ^ijklVkl + 2 ^ijklmnVkl Vmn + ■ ■ ■ - ^k Q^a + ^k g^a + • ■ • ■ ) 

The irreducible components of and have the same form as those of rj" and t/*" , 
and it is i° = where the possible pairs of (a, i) are again evident from Appendix 



A| . Cauchy stresses cr" and cr™ are measured in the distorted state. They follow 
immediately from the Piola-Kirchhoff stresses Because E-^ remains unchanged 
even for a finite deformation, the phasonic case must be treated with some care: 

4 - F,, F,i til , a- = A-i F,, , (6) 

where A = deti^. 

Since our intention is to produce certain strains by means of applied stresses, we 
should rather work with the isothermal Gibbs enthalpy G(t",t"') 

G = Sflf + Sf^lf" + ... . (7) 

S'f'' and Sf'^'^ are independent elastic compliances, and the same invariants as in 
appear, but this time formulated with components of and . One obtains the 
strain-stress relations 

fIC 1 BT'^'' pijabc 

a _____ Q*^^ _ Q<^bc ±b ±c gab ^ k oabc ^ k (o\ 

dtf- ijkl^kl 2 ijklmn^kl ''mn "i • • • ^k Q_^a k Q_i_a ■ ■ ■ ■ \o ) 

The irreducible form hereof is obvious. 



3. The elastic invariants 

First thing of all is to note the transformation behaviour of the strain tensors Q, 
which follows from the transformation behaviour of vectors in i?" and E-^: 

u: (r^ (^r^)^ = ® ^^ w. v^' . (9) 

Index s means symmetrized. The irreducible components of the strain tensors are 



given in Appendix A, and the associated transformation matrices are displayed in 
Table Al. They are deduced from the transformation of the cartesian strain tensors 
'7" = "Hij Vij I where 77°^- is a basis deformation with component ij being 1 and 
all others 0: gij" = vfj 9Vfj = {gV^h Vfj , where (gr/")ij = D^^^{g) D^g) rili , 
{9'n^)ii = Djkjg) ^uid) Vki : ^-iid D^{g) , {g) are the coordinate transformation 
matrices of Table |A1| . 

The characters of symmetrized product representations D — {D ® D)s and 
D = (£) (g) D (g) D)s, respectively, of one and the same representation D, which the 
terms occuring in (H) transform after, are 



X^(.9) = i[x^(5)^ + X^(ff^)] and 

(.9) = i [ X^ (.9)' + 3 x^ (.9^)X^ (5) + 2 x^ (5^)] 



(10) 
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Table 1. Tabulation of symmetrized third-order vector spaces. In the last 
column, all possible third-order expressions consisting of irreducible strains, 



referrin 


g to orthonormal basis sets, are listed. 




Case 


Dimension dim of 


Components with respect 






symmetrized space 


to an orthonormal basis set 




uuu 


56 






uuw 


189 






uww 


270 




U 7^ k) 


www 


165 







Herewith, the foUowing Clebsch-Gordan series for the representations acting in the 
respective symmetrized product spaces of second-order tensors are straightforward 
and given befow in the order uu, uw, ww, uuu, uuw, uww, www |p^ : 

[(r^ ® r^) (g) (r^ ® r^)]^ = e e 3^^ 

(r^ ® r^) (g) (r^ ® r^) = ® 2t^ © 2t^' ® av^ ® 5^^ 

[(r"* ® r^) (g> {V^ ® r^)]^ = 2r^ ® ® r^' ® sr"* ® 5^^ 

[(r^ ® r^) (g> (r^ ® r^) ® (r^ ® r^)], = ® ® r^' ® 4r'^ ® 6^^ (ii) 

[(r^ ® r^) «) (r^ ® r^)], «) (r^ ® r^) = 4ri®8r3®8r3'® lar-*® l7^^ 

(r^ ® r^) ® [(r-* ® r^) ® (r" ® r^)], = rr^® iir^® ur^'® l8^4®25^^ 

[(r4®r5)® (r4®r5)0(r4®r5)], = ®7t^ (b7t^' ®i2t'^®ut'^ . 

We see that we have the fohowing numbers of invariants: 2 (uu), 1 (uw), 2 (ww), 4 
(uuu), 4 (uuw), 7 (uww), 5 (www). These numbers have aheady been calculated in 
an earlier work [T^ . 

The second-order invariants can readily be written as simple scalar products 
^u,i . ^u,i^ ^u,5 . ^u,5^ ^u,5 . j^w,5^ ^w ,i . ^w,i ^w,5 . ^w,5 ygctors Containing 

the irreducible strain components (see Appendix A and jl^). Of course, most of the 

third-order invariants are more complicated. 

In Table |^, the symmetrized third-order vector spaces are specified in more detail. 

Some terms must be weighted with factors to become components for an orthonormal 

basis set and to transform orthogonal among all others. The third-order elastic 

invariants are the components for the basis vectors of the identity representation. 

Basis vectors for an irreducible group representation a are obtained by means of 

the projection operators Pfl = SgsG -^?k* id) (d) with D(g) being the linear 

operators acting on the full (reducible) vector space (see Table |]). These projectors 

Tke^J = Sc.pSjkefi. Here, ef^- 



have the property -Pj'^ef = SapSjksfi. Here, e^- is a basis vector transforming as 



the index j of the irreducible representation /?, and 1 < i < n^, where n^, is the 
multiplicity of /3 in D. To split up the full vector space into subspaces spanned 
by orthogonal irreducible basis vectors, the vectors e°[ S ImP[\, ef^ _L e"^ for 
i ^ j, must be calculated. The other basis vectors are then eg = Pj^ef^, where 
i e {l,...,ng}, j G {2,... da}. The operators D(g) are calculated as orthogonal 
transformation matrices for the components of Table |l]. 

Herewith, the third-order invariants can be found directly from the full spaces of 
Table |l|. However, we picked the third-order invariants by the following procedure: 
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Calculate the second-order irreducible components for the irreducible representations 
r^, and occuring in the first and third series of (11). Then set up all possible 



invariant scalar products with irreducible components of 77" and rj™ and throw away 
occuring linearly dependent invariants. 

Since the components for different irreducible representations do not mix under 
transformation, another very elegant method, leading immediately to the ordering 
with respect to irreducible representations described below, is to expand the triple 
products of 0) completely, for example [(F^ © T^) ^ (F^ ® T^) (g) {T^ © T^)], = 

(r" (g) r^), ® [(r* r"), r^] ® [r^ (r^ o r^),] ® (r^ o r^), , and search 

the invariants for each of the arising triple products of irreducible representations 
separately. 



The third-order invariants are listed in Appendix B. They consist of components 



of as few different irreducible representations as possible (see Table |Bl|). For each 



of the four cases, they are orthonormal, i.e. '^'l^J^Vi,kVj.k = ^ij, where dim are the 
respective numbers in the second column of Table ^ and Vi^k is the coefficient of the 
third-order term k in the invariant i. Furthermore, we have tried on the one hand 
to choose the invariants as short as possible and on the other hand to bring them 
to a suitable for m for compar ing with each other and with the invariants of [|l4[ (see 



Appendix B and Appendix C for details). 



4. Discussion 



From (|11|), there are 20 independent third-order elastic invariants, or elastic constants, 
describing the non-linear elasticity of icosahedral QC. Since there exist four uuw- 
invariants, we have four non-linear phonon-phason couplings. The other third-order 
invariants are unsuitable if one wants to generate phason strains or stresses. Despite of 
this, from all the third-order invariants, the four uuu ones will play the most important 
role because of their influence on the phonon wave propagation. To our knowledge, 
these third-order phonon elastic constants have not been determined for QC so far. 

Generating phason F^-strains from phonon stresses now is possible through the 
invariant /l'""". Due to (||), 

QUUW 

1 , /on To I o ('-i-U\2 o (_lu\2 



- ^ V30 [3 itlY + 3 {tlY - 2 {tlY - 2 {tlY ~ 2 it^Y 



V2 
V3 



OUUW 

-^Vi{2tltl + ^tlt^), 

QUUW 

(12) 

^{u-,ti + Vititt-Vltiti)^ 

OUUlO 

< = - ^ {U^tl - ^tltl + VQtltt). 

Here, a perhaps unexpected factor of 3 is present because of the summation rule in 
(||). From (p^), it is obvious that phason F^'-strains arise from shear stresses ^2 j • ■ • i *6 ■ 
rff is obtained, for example, by applying the stress = t and all other — 0, 772" 
by applying = t, = ^/ij^t, r]^ by applying = t, = ^Jij2t and t]X via 
= = 1/3/2^. The magnitude of an eventually existing phonon F^-stress, 

which is hydrostatic pressure t", has only an indirect effect on the phason F** strains, 
due to the pressure dependence of elastic compliances. Note that the F''-symmetry 
is unlikely to exist without simultaneous phason F'''-symmetry, which is generated by 
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shear stresses even in the hnear, but also in the non-linear regime according to higher 
order compliances (see Table Bl). 

Another possibility to obtain phason F^-strains is, for example, the quartic 
electrostriction [Ea. 



Appendix A. Icosahedral irreducible strains 



The icosahedral irreducible strain components given bel ow are from [ pT[ . They refer to 
the same coordinate systems as in loHI . In Table ^ the icosahedral irreducible 



transformation matrices are given, which the irreducible strains and stresses transform 
after. Other coordinate systems in use are compared to our in some detail in 



V" 



1 



u,5 



W.5 



/ V2 \ 

V J 



/ 2^3 ( 



V 

1 



i%"2 + (T+7te) \ 

^^722 +'733) 



Till 
1 



<3) 



7W1 
7%"2 
Hz 



1 

7& 



72 

%"2 + %"3 

Tr,r2 

V3(r;ri-%1>) \ 

V2(^%"^2-7%"3) 
V V2(r7y^3-7%"i) / 



(A.l) 



Appendix B. Third-order icosahedral invariants 

/«™ ^ _Lyio [(772")' - 3 - 3 %"(r?5")2 - 3 %"(%")2 + 6 ^"(r,^")^ + 3a/3 ^"(t?^')' 

-3V3%"(%-f], 

= iyis^r [(^2)' + + + + 

(B.l) 

= l^/30 [2 (7yJ)2,7r + 2 (%")2< + 2 (Tye")^,;^ - 3 {^^fvT " 3 (%")'< 
+2V5%"%>2"' + ^VEvtVeVs + ^V5V4VbV7 - VTOV3V5V3 
~VT0v3VeV7 + 2VTOV3V4V2 - V30rj^v^V3 + %/30%"776<] , 
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Table Al. Tabulation of the transformation matrices for the icosahedral 
irreducible representations. They are given for two appropriate generating 
elements of Y, i.e. a fivefold rotation C5 and a threefold C3. r — 1 = i. For Yjj, 
the inversion i must be included, i doesn't change the strains 



Cs 



= 1 



= 3 



r^' = 3' 



r5 



T - 1 
1 

1 - T 

— r 
1 



— T T — 1 

-T -1 

1 1 - r 

T — 1 — T 



1 

1 
10 
10 

1 
10 
10 







f 


-1 




-V5 


^1 










/ 


1 










\ 






1 






-1 -3 


-1 

























1 






4 






3 1 


-1 

















1 
















\ 


V5 


1 -1 


-3 












V 








1 




) 






1 


-V^ 


-V& 















-1 














\ 






-1 


-V2 - 


2V2 




V2 












-1 















-v^ 


-V2 


2 







2 




1 

2 




















2 







2V2 





-2 




2 


















2 








v 




V2 


-2 


2 







) 




v 















2 


; 



-3 (773")2< - 4,73"r;4«777- +677"%"%- + 4^2 ^^^^^^ + 4y2,74«%"< 
+A^r^lrilri^ ^V^i^^f^^ - V3 + 2^3 

-2^3%"%"%"'], 



-'3 



+2V3%>,^%»-2V3ry3"%"%-J 

(B.2) 

J™ = + + W)' + W)' + W)'], 

J"- = lV30 [r;^'(%"')' + - 2 7y3"(,7r)2 + 2%"« + 2%"« + 3 

-3%"(C)' - 4r74"C'77"' + 6%"%"'C + 4v^%"%"'< + 4^2 r;«777"'r;- 
+4^/2 ,74"%"'%"' - W)' + V3%"(,7-)2 ~ 2^3%"^%"' 
+2V3%"C%"'], 
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+2x/3«%- - 2^/3 , 



TUWW 



+2y5%"r;2™%"' + 2V5%>2"%- + 2V5n^V3Vr + '^^vlvXv^ 
-%/T0%"« - ^^%"'?^%"' + 2%/1077^%-C + 2%/l0%"%"'??^ 

= 1^30 [ - {rj^f + + - 2 r,6"'(r?7-)2 + 3%"'(%-)2 



(B.3) 



ir'^ = ^vm [2<w)' + 2<w)2 + 2ryrwf - 3 - 3<(cr 

(B.4) 

The invariants comprise strain components of irreducible representations as displayed 

in table Bl. Note that the following pairs of invariants have exactly the same 

structure" /^^^ ^ jwww juuu ^ jwww juuu I j^uww juuw ^ jwww juuw I juww juuw I JUWW 

juww/jwww_ Furthermore, 3 /I*""" is obtained from 7^""" by replacing all products 

VTvfVk by vf-sVfVk + vTVj-sVk + vTvfVk-s- The same is true for 3/^""" and 
j-www _ follows from J™™™ by replacing only the phason F^-components in 

the manner described above. 
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Table Bl. Powers of irreducible representation components in third-order 

invariants. Possible combinations are: a = 1, i = 1; « = 5, i = 2, 3, 4, 5, 6 (a = u); 



o = 4, 


i = 1,2, 


3,4; Q = 


5, i = 5, 


6,7,8,9 (a = wi). 








(a, a) 




{u,5) 


(w,4) 


(w,5) (u,l) 


(u,5) 


{w,4) 


(w,5) 


TUUU 
'l 


3 






JUWW 

^3 


1 




2 


TUUU 

^2 




3 




^4 


1 




2 


TUUU 

^3 




3 




^5 


1 


2 




TUUU 

^4 


1 


2 




^6 


1 


1 


1 


TUUW 
^1 




2 


1 


JUWW 


1 


1 


1 


TUUW 

^2 




2 




-j^ JWWW 




3 




TUUW 

^3 




2 




1 JWWW 

^ ^2 






3 


TUUW 

^4 


1 


1 




1 JWWW 

^ ^3 






3 


TUWUl 
^1 


1 




2 


^4 




2 


1 


JUWW 

-'2 


1 






2 JWWW 




1 


2 



Appendix C. Relationship to other third-order invariants and elastic 
isotropy in phonon space 

tr(r;"){tr[(r7")2]}, = VS/r" + \/5/r" , 

9 12 3d 

det(j7'"') - + -VS/a""""^ - i Vis/™'""' + iVSO/g"'"'"' . 

9 9 9 9 



In (CJ), index s denotes that 77^" must be replaced by ■^{rjfj + rij^i)- Some of our 
third-order invariants are just proportional to 77" times one of the second-order, 
invariant scalar products given in Section ^. Note that tr(r7"') is not invariant. 
Our invariants /I*"" and I^"^^ are also 0(3)-invariants, while J^"" and /^"" must 
be combined to ^^14-^^2"" - ]^\/70/3"" to give a hnearly independent third 0(3)- 
invariant. In the degenerate case of phononic isotropy, our elastic constants (7""" 
can be expressed by sets of classical third-order elastic constants already in use: 

Up to normalisation factors, the phason third-order invariants Xg, X5, J7 and J^' 
of Ishii are transformed into our Hf"""" , 7^'""', I^"""" and /g"'""", respectively by 
the substitutions 77^ ^ -i(%/T07/^ -t- %/67y^) and 7y^ i(-V6ry^ + %/T07/^). This 
is necessary because in other irreducible components are used than in . 
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